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We develop an approach to investigate the non-perturbative dynamics of quantum field theories, in 
which specific vacuum field fluctuations are treated as the low-energy dynamical degrees of freedom, 
while all other vacuum field configurations are explicitly integrated out from the path integral. We 
show how to compute the effective interaction between the vacuum field degrees of freedom both 
perturbatively (using stochastic perturbation theory) and fully non-perturbatively (using lattice 
field theory simulations). The present approach holds to all orders in the couplings and does not 
rely on the semi-classical approximation. 



I. INTRODUCTION 



Lattice field theory represents the only available ab-initio framework, which allows to compute matrix elements of 
a large class of quantum field theories, in a fully non-perturbative way. In particular, due to the continuous advance 
in the development of new machines and new algorithms, lattice calculations for QCD are now beginning to explore 
the chiral regime and are already producing accurate results for a large class of observables. 

On the other hand, lattice simulations do not directly explain the qualitative physical mechanisms which are respon- 
sible for the non-perturbative phenomena. It is therefore important to continue developing alternative approaches, 
which can provide physical pictures and direct insights into the qualitative mechanisms. 

In the specific context of QCD, a large effort has been made in the last decades, in order to identify relevant 
low-energy vacuum gauge field configurations, which are responsible for hadron structure, by driving the breaking of 
chiral symmetry and producing color confinement. For example, instantons have been shown to play an important 
role in the breaking of chiral symmetry [1] and instanton models [2\ have been successfully used to predict physical 
properties of light hadrons (see e.g. [3] and references therein). Similarly, vacuum fields made from monopoles |5j, 
center- vortices [6], merons [7], and, recently, regular gauge instantons |8j have been shown to generate an area law for 
the Wilson loop, hence to produce color confinement. 

Once a set of important low-energy vacuum field configurations has been identified, it is natural to address the 
question whether it is possible to build an effective theory, based on such degrees of freedom. In practice, this corre- 
sponds to deriving an expression for the original generating functional, in which the functional integral is restricted to 
the configurations of the selected family of low-energy vacuum fields, while all other field configurations are integrated 
out and give raise to an effective interaction. 

In the present paper, we take a step in such a direction. The main idea is to use lattice simulations to generate a 
statistically representative ensemble of field configurations. Such configurations are then projected onto the functional 
manifold formed by chosen the family of vacuum field configurations. This procedure is conceptually analog to the 
technique adopted in statistical mechanics to evaluate the free energy, as a function of a set of (order) parameters. 
The result is a new exact expression of the original path integral, given in terms of an integral over the collective 
coordinates of the low-energy vacuum field manifold. 

In order to introduce the formalism and illustrate how the approach works, in this first work we consider the simple 
case of a one-dimensional quantum mechanical particle, interacting with a double-well potential. The choice of such a 
toy-model is motivated by two facts: on the one hand, the relevant non-perturbative vacuum field configurations for 
this system are well known: they are the instantons and anti-instantons, which describe the tunneling between the two 
classical vacua. On the other hand, the simplicity of the model allows us to perform detailed numerical simulations 
and test our method. 

The paper is organized as follows. In section [TT[ we introduce our framework for a generic quantum mechanical 
system. From section [Till we focus on the specific case of the double well problem. In particular, in sections [IV] and 
V} we perform perturbative and non-perturbative calculations of the instanton-antiiinstanton effective interaction. In 
^|we discuss the results of the numerical implementation of this method. 

Then, we shall use path integral Monte Carlo simulations to generate an un-biased ensemble of equilibrium field 
configurations and develop a technique to project such configurations onto the vacuum field manifold. It is important 
to stress the fact that this method does not rely on saddle-point arguments. 
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II. EFFECTIVE INTERACTION FOR THE VACUUM FIELD CONFIGURATIONS 



For sake of simplicity, in this work we shall introduce our formalism for a system consisting of a quantum mechanical 
particle, interacting with an external potential. However, the same method can be applied to quantum field theories 
with arbitrary number of dimensions, as long as they can be formulated on the lattice. 

After performing the Wick rotation to imaginary time, the path integral for the system described by the interaction 
U (x) and corresponding to the boundary conditions 



x[-T/2] = -xi x[T/2] = x f (1) 

is given by 

C [T/2] = 
-T/2] = 

where 



nx[T/2]=x f 

Z[x f ,x z \T] = (x f \e- HT \xi) = / Vx e~ (2) 

Jx\-T/2l=Xi 



T/2 

S[x] = I dt 

-T/2 



i 2 (t) TT , ' 

m — — + U(x) 



(3) 



is the usual Euclidean action. 

Let us consider a generic family of vacuum field configurations (i.e. of paths) x(t; 7) , which depend on a finite 
set of parameters 7 = (71, . . . ,7^) and satisfy the boundary conditions 0. The paths x(t;j) form a differentiate 
manifold M, parametrized by the curvilinear coordinates 71, ... ,7^. 

For every given choice of the parameters 7 it is possible to decompose a generic path x(t) contributing to the path 
integral ^ as a sum of a field configuration 2; (£57), belonging to the manifold M, and of a residual field y(t)\ 

x(t)=x(t;<y)+y(t). (4) 

We shall refer to the field y(t) as to the "fluctuation field". However, in the following we shall never require that the 
vacuum field x(t;j) satisfies the Euclidean classical Eq. of motion (EoM). Hence, both x(t;j) and y(t) represent in 
general quantum vacuum fluctuations. 

Let us now derive a particular representation of the path integral ^ in terms of a set of ordinary integrals over 
the parameters 71, ...,7fc and a functional integral over the fluctuation field, y(r). Since the new representation of the 
path integral contains k additional integrals over dji, . . . , djk, we need to impose k constraints. We choose to enforce 
the k orthogonality conditions 

(y(t) • 3 7i (t,7)) = / ' dty(t) g^(t,7) = 0, i = l,...,fc (5) 

J -T/2 

where the functions g%(t) are defined as 



(6) 



In order to clarify the meaning of the condition ^ we observe that the functions {g 7i (t,7)} i=1 ^ identify the k 

ield 



directions tangent to the manifold M of vacuum fields, in the point of curvilinear coordinates 7 = (71, . . . , 7^). We 
consider only choices of manifold and 7 such that the vectors (|6| define a system of coordinates on the manifold. The 
coordinates (^1, . . . , of a point x(t; 7) are defined as: 

*i[£(t;7)] = (x(t; 7 )^ 7l (t,7)) (7) 

* fc [x(t; 7 )] = (x(t- 7 ).g 7k (t,j)). (8) 

Configurations which lie in a functional neighborhood of the manifold can be projected onto the same system of 
coordinates. The components of such paths x(t) are 

[T/2 

*![x(t)] = 0r(t).<? 7l (£,7))= / dt Xi (t) <? 71 (*,7) (9) 

J -T/2 
[T/2 

* k [x(t)] = (x(t) • g 7k (t^)) = / dt Xi (t) g lk (t^). (10) 

J —T/2 
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Vacuum manifold 



Figure 1: Pictorical representation of the projection of the path x(t) to onto the vacuum field manifold. A path x(t) is 
represented by a point in this picture. The constraints ^ imply that the the fluctuation field y(t) is perpendicular to the plane 
tangent to the manifold in the point of the curvilinear abscissas 7 = 7. 



Hence, the condition ([5| imposes that fluctuation fields y(t) should have vanishing coordinates on the system of 
coordinates defined by the vector {# 7i (£, 7)}i=i, 

Let us now apply a standard technique to implement the k constraints ^ inside the path integral ^ [lOl[TT]. We 
introduce a Faddeev-Popov unity: 

1 = Jd^Jvy fn* (fc) (y(*)-ff7<(*>7)) 7) + -*(*)] *M, (ii) 



which serves as a definition of the functional $>[x]. Note that the integration on y(i) in (11) can be trivially performed 
and one obtains 



(12) 



In particular, we are interested in the value of <$>[x] at the point x(t) = x(t;j) + y(t). If we insert (11) in the original 
path integral ([2|, we obtain, after integration over x: 

/k ry(T/2)=0 / \ 

II d ^ / '' V V n^^-^^7) mt-n)+y(t)]e- S[Htn * +ylt \ (13) 
ti Jy(-T/2)=o V i / 

where the dependence the initial and final points and xy enters implicitly, through the vacuum field x(t;~/). 
The path integral (13) can be formally re- written as 

Z(T; Xi ,x f ) = []] d lt e -* F ^->^\ (14) 

J 1 i 



1 = 1 



where ^(7) is defined as 

r y(T/2)=0 



rvi 1 / \ 

F{ 1 ) = -h\og Vy [W k \y(t)-g 7i (t^)) *[£(*; 7) +l/(t)]e"* 

Jy(-T/2)=0 V i / 



S[£(*; 7 )+3/(t)] 



(15) 



Some comments on what we have done so far are in order. First of all we stress that Z(T; Xi,Xf) can be interpreted 
as the partition function of a system with a finite number of degrees of freedom 71, . . . , 7^. The term ^(71, . . . , 7/e) is 
the analog of the (free) energy in statistical physics an will be referred to as the effective interaction. 
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Let us now address the problem of how to compute F(7i, . . . , 7&), using lattice simulations. Let 
{xi(t) : . . . , XN conf (t)} be a statistically representative ensemble of N con f paths (i.e. obtained by means of lattice 
simulations). The coordinates (<E>i, . . . , of each of such paths are specified by the Eq.s (|9|)-([l0|). Using the 
definition Q and the orthogonality conditions ^ we obtain a set of k non-linear Eq.s for the 71, ... ,7^ variables: 

§ 1 [x(t)] = (x(t).g 11 (t^)) = (£(t;7)^7i(*»7)) = *i(7) (16) 

$ k [x(t)] = (x(t).g 7k (t,7)) = (^; 7 )^ 7fe (t,7)) = ^(7). (17) 

Note that, while the coordinates $1, . . . , on the left-hand-side depend on the path x(t), the functions 
^1(7), . . . , ^(7) on the right-hand-side depend only on the set of collective coordinates (71,..., 7^) and are de- 
termined by the choice of the background field manifold and of the parameter 7. Hence, by solving numerically 
such a system of Eq.s, a value for the curvilinear coordinates can be assigned to each configuration. Repeating this 
procedure for the entire ensemble of lattice configurations xi(t),... ,XN conf (t) one determines the probability density 
7^(71, 7/e), which relates directly to the effective interaction 

F(7i,...,7 fc ) = -i logP(7i,..., 7fc ). (18) 



The effective theory defined by the partition function ( 14 ) allows to perform approximate calculations of the vacuum 
expectation value of arbitrary operators 0(t), 

In fact, if the vacuum manifold contains the physically important vacuum configurations, then 0[x(t)] ~ 0[x(t;j)] 
and 

( ( )> - inLdjte-^,..^ ■ (20) 



We note that, while the partition function ( |14| is independent on the choice of 7 — which specifies the system 
of coordinates on the manifold — the effective interaction F(j) and vacuum expectation values of operators may in 
principle depend on such a parameter. However, such a dependence is generated only by the projection of paths which 
contain very large fluctuations, i.e. lie far from the vacuum manifold. To see this, let us consider the projection of a 
configuration x(t) which lies very close to a point on the vacuum manifold x(t',j f ), i.e. 

WxW-xfatfW-O, (21) 



for some 7'. Then, the projection Eq.s (16)-(17) read: 



(x(t) • g 7l (t, 7)) ~ (x(t; i) • g 7l (£, 7)) = (x(t; 7) • g 7l (£, 7)) (22) 

(x(t)-g 7k (t^))^(x(t-j f )'g 7k (t^)) = (x(t;j).g 7k (t,j)). (23) 

A solution of such set of Eq.s is trivially 7' = 7, for any choice of the projection point 7. If the vacuum field manifold 
captures the physically important configurations, the vacuum expectation values of operators 0(t) will be dominated 
by the configurations in the functional vicinity of the manifold. In the limit in which the relevant configurations are 



only those very close to the manifold, the system of Eq.s (16)-(17) have a unique solution and the expressions (19) 



become independent on the choice of the coordinate system on the manifold, i.e. of the parameter 7. Clearly, this 
condition can be verified by comparing the results obtained projecting onto different points of the manifold. 

In the rest of this paper, we shall illustrate how this method is implemented in practice, in the specific case of the 
one-dimensional quantum double- well problem. 

III. APPLICATION TO THE QUANTUM MECHANICAL DOUBLE WELL PROBLEM 



The discussion made so far has been completely general: Eq.s (14) and (15) hold for an arbitrary choice of the 
potential U(x), of the vacuum field manifold M and of the boundary conditions X{ and Xf. As an illustrative example, 
let us now restrict our attention to the specific system defined by the potential 

U(x) =ma(x 2 -p 2 ) 2 , (24) 
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where m is the mass of the particle. We consider the path integral with periodic boundary conditions (see Eq. [I]) 

Xi = Xf = —(3 ( or equivalently: xi = Xf = +/? ). (25) 

In this specific system, a choice of the effective degrees of freedom is suggested semi-classical arguments. We choose the 
vacuum field manifold to be the one generated by the superposition of N instantons and N antiinstantons. Obviously, 
the choice of the optimal number of pseudo-particles depends on the time interval T. If the barrier is sufficiently high, 
one can fix N from the semi-classical 1 tunneling rate [12] 

2N ~ kT, (26) 

where k is the one-instanton measure 



2(2a) 3 / 2 /? 2 

nh • [ 

The curvilinear coordinates 71 = £1,72 = ?i, • • • , 72iv-i = ^at,72At = Fjv represent the collective coordinates of 
each instanton or antiinstanton, i.e. their positions in the imaginary time axis. In particular, we adopt the so-called 
"sum- arcs ate", which consists in simply adding-up the instanton and antiinstanton fields: 

N 

x S2N t N ) = -P^Yl t f ~ + - > ( 28 ) 

fc=l 

where we have labeled with ti, . . An (?i, . . . , ?at) the centers of the instantons (antiinstantons). The path integral, 
re- written as in Eq.(14) reads 

Z[T- -(3, -f3] = J dt u J di x ... J dt N J di N e"* ntuh,.-,t N ,i N ) ^ 

We note that there are only two choices of such collective coordinates for which the field configuration (28) becomes 
an exact solution of the Euclidean EoM: 

1. When all nearest neighbour instanton-antiinstantons pairs are infinitely separated from each other, i.e. 

\U — ^00, \t i+ i ^00 

2. When all nearest neighbour instanton-antiinstantons pairs are infinitely close to each other, i.e. 

ti tj, 



In the former case, one obtains a dilute instanton gas configuration. In the latter case, all pairs annihilate and the 
field reduces a trivial classical vacuum, i.e. x(t) = —(3. For any other choice of the collective coordinates ti, . . . ,?at, 



the field configuration (28) is not an extremum of the action. 

The relative statistical weight of each configuration in the path integral ^ is provided by the exponential factor 
appearing in Eq. (29), which plays the role of the free energy in the statistical mechanical analogy. Hence, the function 
F(ti, . . . , fjv) expresses the statistical and dynamical correlations between the pseudo-particles, induced by all other 
field configurations in the path integral. For example, in the high barrier limit in which the semi-classical dilute 
instanton gas approximation is justified, one has 



F(ti,...,ttf) 



0(h - h) 0(t 2 - h) • • • 0(t N - t N ) k 



2N 



(30) 



As the height of the barrier is adiabatically reduced, the dilute instanton gas approximation becomes worse and 
worse and eventually breaks down. In this regime, the vacuum fields behave as an interacting liquid and the effective 
interaction F(ti, . . . , iw) deviates from the expression (30) and can be written as 



N 



F(h 



,i N )^YFl A (U-U)+F 2 AI (t 



(31) 



In QCD, where a strict semi-classical analysis cannot be consistently applied, the number of pseudo-particles may be estimated form 
phenomenology or lattice simulations. 
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where FJ> A (F ai ) expresses the two-body instanton-antiinstanton (antiinstanton- instanton) correlations 2 . For very 
low barriers, the average instanton distance becomes smaller than the instanton size, and the pseudoparticles "melt". 
Clearly, in such a regime, instantons and antiinstanton fields no longer represent a good choice of low-energy vacuum 
degrees of freedom. In the remaining of this work, we shall consider systems for which the dilute liquid regime is 
appropriate. 

In order to compute F\ A {t' — t) and F A1 (i! — t) it is convenient to integrate out from (29) all instanton degrees of 
freedom, except those of a single pair of pseudo particles. To this end, we rewrite the path integral as : 



Z[T; -f3, -{3} 



dt\ J dti yj dt 2 dt 2 . . . dt]\f e h J + J dt\ \ dt 2 
\ (J dt' J dt gl A (t' - t) + jdtj dt' gf{t' - t) 



dt\ dti dts • • • dtjsf e h 



(32) 
(33) 



The first term corresponds to the case in which the pseudo-particle of coordinate t is an instanton, while the second of 
coordinate t' is an anti-instanton and g 2 A {t' — t) is the corresponding pair-correlation function. Conversely, the second 
term corresponds to the case in which the pseudo-particle at t is an anti-instanton and that at t' is an instanton. In 
the dilute liquid regime, the functions g IA ^ AI \t' — t) relate directly to F2 A<yAI \t f — t) by 



(34) 



where the proportionality factor is controlled by the density. 

In order to extract the instanton-antiinstanton pair correlation function g^ A we consider the path integral with 
boundary condition Xf = Xi = —(3 and parametrize a generic configuration x(t) using the sum ansatz for an instanton- 
antiiinstanton pair, Eq.(28) 



1 — tanh 



2a/3(t-t 1 ) +tanh V2af3(t-t 2 ) | + y(t) 



(35) 
(36) 



where y(t) is a configuration of boundary conditions y(±T/2) = 0, and t\ and t 2 are the coordinates of the two 
pseudoparticles, in the Euclidean time axis. Conversely, in order to evaluate g AI , one should consider the path 
integral with boundary conditions Xf = Xi = (3 and adopt a vacuum manifold based on the anti-instanton instanton 
pair: 



x(t) = xf(t'MM) + y{t) 



1 — tanh 



2af3(t-t 1 ) 



tanh 



[3(t-t 2 )\}+y(t). 



(37) 
(38) 



Since the two calculations are identical, in the following we shall focus on determining g\ A and the I A suffix will be 
implicitly assumed. 

It is convenient to introduce the relative variables 



X 



l 

2 

t 2 



(*i+t 2 ), 

-h. 



Notice that variable x ^ the "center of mass" of the pair, while £ represents the "relative distance" between the 
instanton and antiinstanton. Notice also that Eq.([33|) implies 

(39) 



F 2 (t 1 ,t 2 )=F 2 (t 2 -t 1 ) = F 2 (0, 



that is to say we expect the effective interaction to be independent from the center of mass of the pair. This is a 
consequence of the time translational invariance of the vacuum. 

We recall that the multi-instanton field configuration and the fluctuation field have to fulfill the orthogonality 
conditions ([5]), which is enforced in a specific point 7 = 7 of the manifold. The basis vector of the tangent space of 



Eq. |3l| can be generalized to include higher-order (e.g. three-body, four-body, etc..) correlations. 
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the manifold defined by the sum ansatz (28) are, for an arbitrary point 7 = (£1,^2) 
g tl (t'MM) = d tl xs 2 (t]t 1 ,t 2 ) 

gt 2 (t;ii,i 2 ) = d t2 x S2 (t;t u t 2 ) 

t\=t\ ,t 2 =t 2 



t\=t\ ,t 2 =t 2 



2o^ 2 sech 2 



2o^ 2 sech 2 



2af3(t - i 2 ) 



(40) 
(41) 



Equivalently, in terms of the x an d £ coordinates, the basis vectors of the tangent space in the generic point 7 = (£, x) 
read 



1 



1 



d x %s 2 



= V^/3 2 < sech : 



2a/3 t - x 



x=x,€=€ 



sech 



2a/3 [t - x + 



sech 



2a/3 (t-x 



sech 



2a/3 U - x + I 



Hence, without loss of generality, in the following we shall consider 



z(t;x,0 

with the conditions 



/3 < 1 - tanh 



2a/3 [t-x + 



tanh 

= 0, 
= 0. 



2a (* " X - |) } 



■!/(*) 



(42) 

(43) 
(44) 

(45) 
(46) 



Although our ultimate goal is to evaluate an d ^/^(O m a fr^Y non-perturbative way, it is instructive 

to discuss first a perturbative analysis, which yields information about the contribution to the quantum effective 
interactions in the short instanton-antiinstanton distance limit. Such a calculation isl be presented in the next 
section, while the fully non-perturbative calculation is reported in section |V| 

IV. PERTURBATIVE CALCULATION 



Perturbation theory deals with small quantum fluctuations around a classical vacuum. In particular, a calculation 
of i^2 /A (0 an d ^2 /A (0 a ^ small £ requires to assign to each point in the vicinity of the trivial vacuum 



x S2 



(47) 



a point on the intanton-antiinstanton functional manifold. Since quantum fluctuations can be arbitrarily small, the 
orthogonality conditions (45) and (46) have to be imposed at a point which is arbitrarily close to the same classical 



vacuum. In principle, the most natural choice would be impose the orthogonality conditions at the classical vacuum. 
However, problems arise due to the fact that it is not possible to define the tangent space in such a point, since 



g x (t;x,0) = 0, 



(48) 



for all x- To overcome this difficulty, in the following we use the stochastic quantization formalism to construct a 
rigorous approach in which the tangent space which is defined at a point which is arbitrarily close to classical point, 
but does not coincide with it. 

Let us begin by briefly reviewing Pairsi and Wu quantization technique [14J. The starting point is to allow the field 
configuration x(t) to depend on an additional parameter, the so-called stochastic "time" r. The dynamics of the field 
in such an additional dimension is postulated to obey a Langevin equation: 



(49) 
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where k is an arbitrary diffusion coefficient and rj(t, r) Gaussian distributed stochastic field 

P[rj\ (x expj-^ J dtj^ drr/ 2 (t,r)J 
which obeys the fluctuation-dissipation relationship 

( V (t, r)rj(t', r')> = J Vf] v(t, T) V (t', r')P[ V ] = 2 k 6(t' - t) 5(r' - r). 



(50) 



(51) 



For any value of the stochastic time r, the probability to for the field to assume a given configuration x(t,r) is 
described by a (functional) probability distribution V[x](r), which is a solution of the Fokker-Planck Eq. associated 
to the Langevin Eq. (49): 



d ^ n , S 2 ^ . 1 S / . _ SS[x] 



(52) 



A general property of the Fokker Planck Eq. is that its solutions converge to the static, "Boltzmann" weight, in 
the long time limit: 



(t— >oo) 



JVx(t)exp{-j i S[x(t)}} 



exp 



{-\s[x(t)]}, 



(53) 



regardless of the initial condition, x(t,r = 0) and of the value of the diffusion coefficient k. Hence, the Langevin 
Eq. (49) generates configurations which, at equilibrium, are distributed according to the statistical weight appearing 
in the Euclidean quantum path integral. Such configurations can be used to compute quantum mechanical Green's 
functions. 

In stochastic perturbation theory, a generic path x(t,r) obeying Langevin Eq. (49) with boundary conditions 
is written as a power series in e = \fh\ 



x(t,r) = ^2e l Xi(t,r), 



(54) 



xo(t,r) is the classical content of the path, while all other terms represent quantum corrections. In the double- well 
problem, the classical solution with boundary conditions (25) is xo(t,r) = —j3. 

By inserting the expansion (54) into the Langevin Eq. (49) and matching the Left-Hand-Side (LHS) and Right- 
Hand-Side (RHS), order by order in e, one generates a tower of coupled stochastic differential Eq.s, for the components 
Xi(t, r), which appear in Eq. (54): 



0(e°) : 


x' (t,r) 


= km (d 2 - 


-4a [x 2 (t,r)-p 2 ] ) 


x (t,r) 






Oie 1 ) : 


x[(t,r) 


= km (d 2 - 


-4a [3x 2 (t,r) - p 2 } 


) zi(*,t)H 






0(e 2 ) : 


x 2 (t,r) 


= km [(d 2 


-4a [3x 2 (t,r) - p 2 \ 


) x 2 (t,r) 


— 12a xo(t, r] 


)x 2 (t,r)] 


0(e*) : 


x 3 (t,r) 


= km[ [d 2 


-4a [3x 2 (t,r) - p 2 \ 


) x 3 (t,r) 


- 4a x\(t,r) 


— 24axo(t, r)xi(t, r)x 2 {t^ r 



(55) 
(56) 
(57) 
(58) 



In practice, the perturbative expansion is truncated and one solves a finite set of stochastic differential Eq.s, starting 
from a given initial condition. For example, truncating the expansion to order e 2 and choosing the initial condition 



z (t,T = 0) = -A 

X z (t,T = 0) = 0, (2 = 1,2,...), 

which corresponds to the classical vacuum state, we find 
x (t,r) = -p 



Xi 



(t,r) = r -e-** f 

J-oo 27r JO 



dr 



-fcm(8a/3 2 ^)(TV)~^ T ') 



duo 



x 2 (t,r) = Ylapkm I —e~™ 1 \ dr 

27r Jo 



dt 



X ((t,T). 



(59) 
(60) 

(61) 
(62) 

(63) 
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The corresponding perturbative solution is 

x(t,r) = x (t,r) + exi(t,r) + € 2 x 2 (t,r). 



(64) 



It is important to stress that only the asymptotic equilibrium solution x(t,r = oo) enters in the evaluation of physical 
observables. Such equilibrium solutions are independent on the choice of the initial condition of the perturbative 
stochastic equations (55)-(58). 



Let us now show how the stochastic perturbation theory technique can be used to gain information about the 



Fl A (0 an d ^2* 7 (0 distributions. To this end, we begin by decomposing the field as in Eq. (|4|, 

x(t) = xs 2 (t;t u t 2 ) + y(t). 



(65) 



Next we need to promote the manifold field x(£;x,0 and fluctuation field y(t) to dynamical variables, under the 
stochastic time evolution. There is some freedom associated to the definition of such a stochastic dynamics. For 
example, a possible choice may be one in which the r dependence enters entirely through the fluctuation field y(t,r), 
while the smooth vacuum field xs 2 is assumed to be static, under stochastic evolution, i.e. xs 2 (t 1 r) = xs 2 (t). Instead, 
a crucial point of the present approach is to make a different choice and allow both the fluctuation field and the smooth 
vacuum field to vary with the stochastic time r. This is done in practice by promoting the curvilinear coordinates 
£ and x to dynamical stochastic degrees of freedom [15], i.e. £ — > £(r) and x ~^ x( r )- Consequently, at a generic 
stochastic instant r, the quantum field x(t,r) reads: 



x(t,r) = x S2 t;x(r) 



^(r),x(r) + ^(r))+y(t,r). 



(66) 



Let us now construct a perturbative solution of the Langevin Eq. (49), based on the decomposition (66). We 



recall that the multi-instanton field is not a classical solution of the EoM, except in the points where £ = 0. As a 
consequence, quantum corrections will appear not only in the fluctuation field, but also in the background field. To 
account for this fact, we expand y(t,r), %(t) and £(r) as power series in e = ^/h: 



y(t,r) 

X(r) 
C(r) 



oo 
i=0 

CO 



(67) 
(68) 
(69) 



i=0 



Let us now define the tangent space in a generic point £, x of the manifold. It is possible to show that the orthogonality 
conditions (45) and (46) hold order-by-order in perturbation theory and at any stochastic time i.e.: 



Vi(t,T) -g x (t;x,0j = 0; 
^(^).^;x,0) =0, Vi,Vr. 



(70) 
(71) 



From Eq.s (54), (67), (68) and (69) it is immediate to obtain an expression for each of the Xi(t,r) components in 
Eq. d54b: 



1 d i 
i\ de i 



•E S<2 



n=0 



m=0 



1=1 



(72) 
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For example, the first orders are 
x (t,r) 

x 2 (t,r) 



xs 2 U;xo ^,Xo -' 



2 ,/vu 2 

XiW #x(^;Xo,£o(V)) + £i(t) ^(t;Xo,Co(r)) + 2/i(*,t) 
X2(r) # x (£;Xo,£oO)) +6W ^(^Xo^oW) + 



4 (xi(r) + ^ ) ^(t,Xo,eoW) -XiW 6(r) Xo, CoM) + y 2 (t, r) 



(73) 
(74) 

(75) 



where we have used the fact that xo is independent on r. The terms on the LHS of Eq.s (1731- (75) coincide with the 



perturbative solution results (61)- (63). On the other hand, the terms on the RHS represent the decomposition of the 
same functions in terms of the low-energy vacuum field configurations and of the corresponding fluctuation fields. 

In order to make contact with the effective interaction, we need to introduce the tangent space which enters the 
projection Eq.s (70) and (71). At this point, we need face the above mentioned problem that the tangent space at 



the classical vacuum —ft is not defined. To overcome this problem, we let the tangent space vary with the stochastic 
time in such a way that the point C, x asymptotically approaches the classical vacuum, but does not coincide with it 
at any finite r. In practice, we promote £ to a stochastic variable and we impose 



0. 



(76) 



In particular, we choose £(r) = Co( r ), since Co( r — > oo) — > 0. 

Using such a decomposition, we are now in a condition to analytically compute arbitrary moments of the equilibrium 
distribution for x an d C, i.e. (C fc ) and (x k )- By projecting and inverting Eq.s (74) and (75), we obtain the following 
expression for the collectives coordinates up to 0{K) 



Xi (r) 

X2(r) 
6(r) 



(xift,r) •g x (t',xo,to(T))) 

Xo, fo(r) • g x (t; xo, foM) ' 
(xi(t,r) •^(t;xo ? go(r))) 
(* ; Xo , Co (t) • ^ (t; Xo , Co (r ) ) ' 
(x 2 ft,r) -p x (t;xo^o(T))) + Xi(t)£i(t) (gjftj Xo, CoCO) ; 9 X (^ Xo, £o(t))) 

(#x Xo , Co (r ) ) • g x (t ; xo , Co (r ) ) ) 
MVr) •^ft;xo,Co(T))) + | (Xi(^) + |CiQj) fefe Xo, CoQj) ; g^j Xo, go(r))) 

Xo , Co (r ) • ^ (t ; xo , Co (r ) ) 



Using the fluctuation-dissipation relationships (51), and the fact that Co(^) is independent from rj(t,r) we find 



<Xi(r)> 


= 


<X2(r)> 


= 


(6(r)> 


= 


(6(r)> 


_ ( 



(&(*; Xo,£o(t) •#(i;Xo,£o(T)) 



T — > OO 

9 1 



32 ma/3 4 

Hence, we have obtained a closed analytical expression for the first moments: 

(x> = (Xo + £Xi +£ 2 X2> - " 1 ruh2 
(0 = <£o + ££i+£ 2 6> = 



= Xo + 0(^ 2 ) 
9 ft 
32 maf3 4 



(77) 
(78) 
(79) 

(80) 
(81) 



(82) 
(83) 
(84) 



(85) 

(86) 
(87) 
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Now, in order to compute the second moments, we observe that the general expression up to order 0(e 2 ) is 

X 2 (t) = xl(T) + 2ex Xi(T)+e 2 [2xoX2(T) + xi(T)] + ... (88) 

C 2 (r) = eo(r) + 2e^(T^ 1 (T) + e 2 [2^(T)^(T)+e i (T)} + ... (89) 

which immediately gives 

(X 2 ) = oo; (90) 

Some comments on these results are in order. The distribution of the instanton-antiinstanton distance £ is not 
symmetric around the origin, since (£) 7^ 0. Such a symmetry breaking comes from fluctuations which explore the 
non-harmonic region of the potential function U(x). Since the potential on the left of the equilibrium configuration 
raises more steeply than that on the right, quantum paths in the direction of the barrier are statistically favored. The 
divergence (x 2 ) = 00 emerges because the distribution of collective coordinates is independent on x? as consequence of 
the time-translational invariance of the system. In the language of stochastic quantization, this implies that the center 
of mass of the instanton-antiinstanton pair performs Brownian motion in stochastic time and (x 2 ) oc r, according to 
Einstein relationship. 

We emphasize once again that in this calculation we have never requested that the multi-instanton configurations 
should be approximate solutions of the classical EoM. The only request is that the configuration corresponding to 
the classical vacuum must belong to the manifold parametrized by the relevant low-energy degrees of freedom. We 
also stress the fact that there is no contribution to the effective interaction, at the classical level: 1^(0 is an entirely 
quantum effect. 



V. NON-PERTURBATIVE CALCULATION 



Let us now take the main step of the present work and perform a fully non-perturbative calculation of F 2 (£) which 
describes the correlations between consecutive tunneling events. 

Let r = 00), . . . , x l (t, r = 00)} be an ensemble of I equilibrium field configurations, which were obtained 

non-perturbatively, for example by integrating numerically directly the Langevin Eq. (49), or by means of a lattice 
Monte Carlo simulation. The pair correlation function gi A {C) can be extracted by projecting the set of equilibrium 
configurations onto the low-energy vacuum field manifold spanned by an instanton-antiinstanton pair. To this end, 
we define the functionals of the field configuration x(t, r) 



$ x [x T ] 
^[x T ] 



(x(t,r),g x (t',x,€)) , 
(x(t,r),^(t;x,0) , 



(92) 
(93) 



which represents the projection of an arbitrary field configuration onto the tangent space, at the point (x,£). We also 
introduce the functions of the collective coordinate x an d £ 



*x(x,0 



XS2 \ t, X 
XS2 [t,X 



\^x+\i 



(94) 
(95) 



which represents the projection of a generic point of the instanton-antiinstanton field manifold onto the same tangent 
space. In the specific case of the double-well potential one has 



*x(x,0 



C 



x 



x ■ 



-2(3' 



c 



x + 



x ■ 



x- 



x + 



x + 



x ■ 



(96) 



(97) 
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Test of the calculation of F 2 ( 

Dilute instanton gas 

4 | , 1 , 1 , 1 , 1 , 1 , 1 , 1 , 1 , 1 , p 

Calculated 

Expected theoretical result (dilute gas limit) 




Figure 2: Non-perturbative calculation of the effective interaction i*2(£) for a dilute instanton gas with a = 7,m=l,f3=l 
in a volume T = 200. The points are the results obtained from projecting 1000 configurations, while the dashed line is the 
expected theoretical results for £ much larger than the instanton size (which is 0.26, in these units). 



where 



C[X] = 2p 2 jv^XsinlT 



2apX 



coth 



2apX }. 



(98) 



By setting Eq.s ( |92| and (93) to be equal to * x and *j respectively, we obtain a complete system of equations for 
the variables \ anc 



* x [ar T 



*x[x(r),C(r)] 
$ € [x T ] = * C [x(t),£(t)]. 



(99) 



Such a system has a unique solution for any choice of the projection point (x,£), with £ ^ 0. Hence, it is possible 
to assign a value of x an d £ to every non-perturbatively generated configuration x(t,r). 

Repeating such a projection for the entire ensemble of equilibrium configurations, one obtains an histogram which 
by construction is proportional to the pair correlation function #2^(0- The effective potential i^(0 is immediately 
extracted from: 



H A iO = "ft Iog[ gi A (Z)] + const.. 



(100) 

Clearly, the calculation of F^ 1 (0 would be completely analog. In practice, such a calculation is not necessary, since 
the function F 2 A7 (£) can be inferred directly by symmetry arguments: 



F 2 AI (0=Fl A (-H). 



(101) 



Once the effective interaction has been determined, one can evaluate the instanton density of the liquid by minimizing 
the free-energy of the ensemble. In the next section, we present the results of some numerical investigations, in order 
to illustrate the method and assess the accuracy of the determination of the instanton-antiiinstanton interaction. 



VI. NUMERICAL STUDIES 



In order to show that our method yields the correct result, let us first consider a model for which the effective 
interaction can be evaluated analytically. The two-body part of the effective interaction for a dilute instanton gas 



can be easily computed from Eq. (30), by integrating out all the collective coordinates, except for those of a single 



instanton-antiinstanton pair. The result is 

e -**i IA K> = const, x e -(«T-i)iog(T-o (102) 
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Figure 3: Non-perturbative calculation of the effective interaction F(£) for a = 1 (left panel) and a = 7. 



This result holds for high barriers and distances £ much larger than the instanton size. Notice that, in the thermo- 
dynamic limit — TV, T — > oo and N/T = k fixed — the effective interaction ^2^(0 should scale linearly, with a slope 
controlled by the instanton rate n. 



We now address the question if our projection technique is able to reconstruct the effective interaction in Eq.(102). 



To this end, we have generated an ensemble of 1000 dilute gas configurations, by randomly sampling the positions of 
instantons and anti-instantons, in a box of size T = 200 for a well with a = 7, m = 1, /? = 1. In Fig.[2j we compare the 
expected theoretical curve (dashed line) with the result of our numerical calculation (points). We see that, as soon 
as the distance £ becomes larger than few instanton sizes — which is 0.26 in this units — the numerical results agree 
with the expected curve. A linear fit of the data for £ > 1 yields a slope of 0.32 ± 1, in excellent agreement with the 
exact theoretical result, which is 0.31. Hence, we conclude that our projection method is indeed able to quantitively 
reconstruct the structure of the exact distribution used to generate the ensemble of configurations. 

Let us now discuss for completeness the structure of the effective interaction, for our original quantum double-well 
system. At this level, we no longer consider the semi-classical dilute gas model. Instead, we account for quantum 
fluctuations to all orders. As the barrier becomes higher and higher, performing a sampling of multiple barrier-crossing 
paths contributing to the functional integral with dynamical algorithms such as Molecular Dynamics of Monte Carlo 
becomes highly inefficient 3 , and computationally expensive. To cope with this problem, we have evaluated the 
instanton-antiinstanton interactions using the importance sampling approach described in appendix [A] 

Fig. [3] shows the results of such a non-perturbative calculation for a well with a = 1 (low barrier) and a = 7 (high 
barrier Some comments on these results are in order. First of all we note that the minimum of the effective interaction 

A (0 18 located at positive values of £, in qualitative agreement with our perturbative calculation. The range of £ in 
which the effective interaction -F/^CO 18 n °t flat corresponds to close, largely overlapping instanton-antiinstanton pair 
configurations. When the distance becomes of the order of twice the instanton size, the effective interaction starts to 
raise and eventually reaches the dilute gas limit. On the other hand, for low barriers, the instanton density is large 
and the attraction and repulsion generated by -F/^HO become important. In such a regime, the vacuum behaves like 
a one-dimensional liquid, rather than as an ideal gas. We note that this is precisely the physical picture underlying 
the instanton liquid model of the vacuum [13J. 



VII. CONCLUSIONS 



In this paper, we have presented an approach which allows to express quantum mechanical path integrals, in terms 
of few ordinary integrals over a set of low-energy variables, which parametrize the manifold of the relevant vacuum 
field configurations. We have developed a rigorous technique to extract the effective interaction, a simple quantum 
mechanical problem, in which the low-energy degrees of freedom are multi-instanton configurations. 

We have assessed the accuracy of our method by correctly reconstructing the effective interaction used to generate 
an ensemble of sinthetic configurations. We have also performed both perturbative and non-perturbative calculation 



3 This difficulty is not present in QCD, where it has been shown that typical lattice configurations contain indeed many instantons and 
anti-instantons. 
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of the quantum effective interaction of an instanton-antiinstanton pair. In both cases, we have found that the effect 
of quantum fluctuations is to shift the location of the minimum of the effective interaction i 7 ^ A (0 t° the right. 

We stress the fact that, although the present discussion has focused on an instanton liquid picture of the vacuum, 
our projection method does not rely at all on the semi-classical approximation. The semi-classical approximation 
has been used only as a guidance to find good vacuum effective degrees of freedom, for our toy model. Hence, the 
method can in principle be generalized to build effective theories for the vacuum, based on different types of field 
configurations. This observation may become important in QCD, where fields other than singular gauge instantons 
are needed, in order to account for confinement. 

If the vacuum fields selected are the ones driving the system's non-perturbative dynamics, then one expects that 
the contribution coming to the fluctuations around them to the field operators appearing in the Green's functions will 
be small. In this case, the calculations of observables in the effective theory can be performed very efficiently, because 
they involve only few ordinary integrals over the set of curvilinear coordinates. 

Most importantly, once a specific choice of the vacuum manifold has been identified, the corresponding effective 
theory yields parameter-free predictions. Hence, the present framework can be used to assess the importance of 
different families of vacuum fields, by directly comparing the results of the corresponding effective theory with the 
experimental data. 

The extension of the present formalism to QCD is in progress and will be presented in our upcoming papers. 

Appendix A: ALGORITHM USED IN THE EVALUATION OF F 2 IA (0 

We are interested in computing numerically the integral 

gi A (0 = M J z>yex P j-| dt£[i Sa (t; X -^,x+^0 + y(*)]}*(tf-flx(x,^)«(y-flc(x,€) 



9x(x>£) -9 x {x,^j (#?(x,£) -9i{x,^j ~ (#x(x,£) (#«(x,£) '5x(x>0 



(Al) 

where the term inside the square brackets is the explicit representation of the Jacobian factor <f>[y]. 

The meta-stability of the double well system makes it rather computationally challenging to generate a statistically 
significative ensemble of field configurations, using algorithms based on Metropolis or by Langevin dynamics. The 
main problem is that, for such a meta-stable system, ergodicity is reached only in an exponentially large computational 
time. 

The problem has no easy solution within a dynamical Monte Carlo approach. However, because of the low di- 
mensionality of our system, simpler importance sampling technique are available and efficient 4 . Since the system is 
time-translationally invariant, without loss of generality we can set x — 0> X — — 1 - e - we can remove completely the 
dependance from the center of mass — , and set £ = 0. Then, the resulting expression for the pair correlation function 
can be re-written as 



9l A {0 =M Vy S(y(t) • gz(t; 0)) <%] P[y(t)} g— , (A2) 

where P[y(t)] is a probability distribution to be defined below. We stress that now the integral can be rest ricte d to 
the small region in which the projection function is not exponentially small. The discretized version of Eq. (A2) is 

gl A (0=^ J l[ay(t k )5 \^y(t k )gz(t k ,0)\ <%] P[y] j- , (A3) 



k=l \k=l 



where N is the number of points in the lattice and Si at is the discretised version of S. 
For P[y] we choose: 



P[y] <* exp [ -— ^ £>(**+i) - Vitk)) 2 ] (A4) 
V k=o J 



4 Note that this problem is unrelated to our projection approach, which is a prescription about the measurement of an effective action, 
once a significant sample of configurations has been provided in some way. 
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with the constraint y(to) = ?/(£jv+i) = 0. We eliminate the delta function by setting the last coordinate y(tjsr) equal 
to 

N-l 

y(tN) = J2 y(ti)gz(U)/gz(t N ). (A5) 

i=l 

Notice that, in this way, the orthogonality condition is satisfied configuration by configurations. The statistical weight 
of resulting each paths was evaluated from 

= exp(-Sutfe(* fc ,Q + w]) j (A6) 
P[Vi] 

Up to an overall multiplicative factor, g IA (x) can be extracted from 



9l Al 



l (0 = const, x ^^(O (A7) 

i 

By taking the logarithm, one obtains Frf up to an overall additive constant. 
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